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Abstract 

■^ ' We investigate the convergence properties of a perturbation method proposed some time ago and 
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<^ ' reveal some of it most interesting features. Anharmonic oscillators in the strong-coupling limit 
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prove to be appropriate illustrative examples and benchmark. 
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I. INTRODUCTION 

Some time ago Bessis and Bessis (BB from now on) 1^ proposed a new perturbation 
approach in quantum mechanics based on the apphcation of a factorization method to a 
Riccati equation derived from the Schrodinger one. They obtained reasonable results for 
some of the energies of the quartic anharmonic oscillator by means of perturbation series of 
order fourth and sixth, without resorting to a resummation method. In spite of this success, 
BB's method has passed unnoticed as far as we know. 

The purpose of this paper is to investigate BB's perturbation method in more detail. 
In Section O we write it in a quite general way and derive other approaches as particular 
cases. In Section lllll we carry out perturbation calculations of sufficiently large order and 
try to find out numerical evidence of convergence. One-dimensional anharmonic oscillators 
prove to be a suitable benchmark for present numerical tests. For simplicity we restrict 
to ground states and choose straightforward logarithmic perturbation theory instead of the 
factorization method proposed by BB PJ. Finally, in Section HVl we discuss the results and 
draw some conclusions. 

II. THE METHOD 

In standard Rayleigh-Schrodinger perturbation theory we try to solve the eigenvalue 
equation 

H^ = E^, H = Ho + XH' (1) 

by expanding the energy E and eigenfunction \E' in a Taylor series about A = 0. This method 
is practical provided that we can solve the eigenvalue equation for A = 0. 

In some cases it is more convenient to construct a parameter-dependent Hamiltonian 
operator H{l3) that one can expand in a Taylor series about (3 = 

H{(3) = Y: H,P'^ (2) 

3=0 

in such a way that we can solve the eigenvalue equation for H{0) = Hq. In this case we 
expand the eigenf unctions \E'(/5) and eigenvalues E{[3) in Taylor series: 

oo oo 

^m = Y.'^,(3\E{(3) = Y.E,(3^ (3) 

i=o i=o 



There are many practical examples of application of this alternative approach p]. In par- 
ticular, BB yy suggested the following form of H{(3): 

W{f3) = T.W,(3=- (4) 

j=0 

Comparing equations ^ and (jlj) we conclude that 

Ho = H- Wo, 

H, = W,.,-W,,j>0. (5) 

In principle there is enormous flexibility in the choice of the operator coefficients Wj as we 
show below by derivation of two known particular cases. 

If we restrict the expansion (jlj) to just one term W{P) = W{0) = Wo then H{P) = 
H — Wo + PWq. Choosing Wo = H — Ho{oi) = H'{a), where Ho{a) is a parameter- 
dependent Hamiltonian operator with known eigenvalues and eigenf unctions, we obtain the 
method proposed by Killingbeck ^3:] some time ago. The main strategy behind this approach 
is to choose an appropriate value of the adjustable parameter a leading to a renormalized 
perturbation series with the best convergence properties l2, IS]- 

If we consider two terms of the form Wq = H — Ho{<y) and Wi = XH', then we derive 
the Hamiltonian operator H{i3) = Ho{a) + l3[Ho — Ho{a)] + j3^XH' that Killingbeck et al 
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have recently found to be even more convenient for the treatment of some perturbation 
problems. Those approaches are practical if we can solve the eigenvalue equation for (3 = 
and all relevant values of a. 

We should mention that it was not the aim of BB to obtain renormalized series with an 

adjustable parameter but to choose the operator coefficients Wj in such a way that they 

could solve the perturbation equations 

j 
(Ho - Eo) ^, =J2{E^-W^ + W^,) ^,_, (6) 

4=1 

in exact algebraic form jjj. 

For simplicity in this paper we concentrate on a one-dimensional eigenvalue equation of 
the form 

^"(x) = [U{x) - E]^{x), U{x) = V{x) + iiilil. (7) 



If ^(0) = \l/(oo) = and / = 0, 1,2, . . . is the angular-momentum quantum number, this 
equation apphes to central-field models. If ^(— oo) = ^^(cxd) = and / = —1, we have a 
one-dimensional model. In particular, when V{x) = V{—x) then / = —1, and / = 0, select 
the spaces of even and odd solutions, respectively. In any case the regular solution to the 
eigenvalue equation (jTj) behaves asymptotically as x^+^ at origin. 

In order to calculate perturbation corrections of sufficiently large order by means of BB's 
method we define 

fix) = --^,s = l + l (8) 

that satisfies the Riccati equation 

f' + -f-f + V-E = 0. (9) 

X 

The corresponding equation for the Hamiltonian H[f3) in equation (^ reads 

f'+'^f-f + V-E + {P-l)W = 0, (10) 

if we restrict to the case that W{l3) = W{P, x) depends only on j3 and the coordinate. The 
coefficients of the expansion 

oo 

f = T.fjP' (11) 

i=o 
satisfy the perturbation equations 

/; + -fj - E M-k + vs,o - E, + iy,_i - w, = 0. (12) 

^ fe=0 



III. RESULTS 

Simple one-dimensional anharmonic oscillators V{x) = x^ + Ax^^, K = 2,3, . . ., are a 
suitable demanding benchmark for testing new perturbation approaches. We easily increase 
the degree of difficulty by increasing the values of the coupling parameter A and anhar- 
monicity exponent K. BB applied their method to the first four energy levels of the model 
with K = 2 and several values of A, restricting their calculation to perturbation theory of 
order four and six. Here we consider the strong-coupling limit (A — ^ oo) of the oscillators 
mentioned above: 

V{x) = x^''. (13) 



Notice that if the perturbation series gives acceptable results for this case, then it will cer- 
tainly be suitable for all < A < oo. Moreover, the perturbation corrections for these models 
are simpler enabling us to proceed to higher orders with less computational requirement. 

In order to make present discussion clearer we first illustrate the main ideas of the method 
with the pure quartic oscillator K = 2. We try polynomial solutions of the form 

/,(X) = J2 Cj,2m+ix'"'^\ J =0,1,... (14) 

m=0 

in the perturbation equations (jT^ for the ground state (s = 0). Substitution of /o(a;) into 
the perturbation equation of order zero leads to 

- cl^x^ + (1 - 2co,iCo,3)x^ + (3co,3 - cli)x^ + Co,i -Eo-Wo = 0. (15) 

In order to have a solution with Co,3 7^ we choose Wq = —cl^x^; then Co,3 = l/(2co,i), and 
3co,3 — c^i = becomes a cubic equation with two complex and one real root. If we select 
the later we finally have 

Wo = -12-2/3a;6^ 

12^/3 
Eo = ^- ^ 1.1447. (16) 



We expect the resulting unperturbed wavefunction 
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^0 oc exp ^-^x^ - ^-^j (17) 

to be an improvement on the harmonic-oscillator one in standard Rayleigh-Schrodinger 
perturbation theory 2, S 1^ • The zeroth-order energy (fT^ is reasonably close to the exact 
value shown in Table IJ which was obtained with the Riccati-Pade method p . 
At first order we have 



122/3 ^ ^,„/ 12l/3ci3 1 \ 

^c,,.«-12V3U + ^^ + - X 
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+ f 5ci,5 - 12V3ci,3 - ^^^] x' + (3ci,3 - 12^/^ci,0 ^' 



6 
+ ci,i -Ei-Wi = 0. (18) 

We easily solve this equation if Wi = —12^/'^ 01^5X^/6; the result is 

^' 112 ^3 56' 



122/3x8 

Wi = , 

112 ' 

El = 12^/^ ^ -0.1022. (19) 

112 ^ ^ 

The energy corrected through first order is somewhat closer to the exact value: Eq + Ei ^ 
1.0425. 

The systematic calculation of perturbation corrections of larger order offers no difficulty 
if we resort to a computer algebra system like Maple p. Since we are unable to prove 
rigorously whether the perturbation series converges, we resort to numerical investigation. 
Figure HI shows that \og\Ej/EQ\ first decreases rapidly as j increases, but then it increases 
slowly suggesting that the series does not converge. If we assume that the error on the 
energy estimated by the partial sum 

M 

eM = E^. (20) 

i=o 

is proportional to the first neglected term \E — E^^^^\ ^ ji^M+il, then it is reasonable to 
truncate the perturbation series so that I-Eaz+iI is as small as possible |7|. In this case we 
find that -E26 = —0.3897686104 x 10~^ is the energy coefficient with the smallest absolute 
value so that our best estimate is E^"^^^ = 1.06036215 (compare with the exact value in 
Table HJ. 

We proceed exactly in the same way for the pure sextic oscillator K = 3. Figure Q shows 
values of log \Ej/Eo\ that clearly suggest poorer convergence properties than in the preceding 
example. The energy coefficient with the smallest absolute value is -E15 = 0.2759118288 x 
10^^, and our best estimate E^^'^'^ = 1.14470 is reasonably close to the exact eigenvalue in 
Table HI In principle, it_js not surprising that perturbation theory yields poorer results for 
K = 3 than ior K = 2 

For the pure octic anharmonic oscillator f^ = 4 we look for polynomial solutions of the 

form 

i+3 

/,(X) = J2 C,,2m+ix''"'^\ J = 0, 1, ... . (21) 

m=0 

In this case we have calculated less perturbation coefficients because they require more 
computer memory and time. Surprisingly, the values of log \Ej/Eq\ in Figure [T] suggest that 
the perturbation series for K = 4 exhibits better convergence properties than the one for 
K = 3 just discussed. The energy coefficient with the smallest absolute value (among those 



we managed to calculate) is: E12 = —0.5205493999 x 10~^ so that our best estimate is 
^[11] = 1.225822 which is quite close to the exact one in Table HI 

The surprising fact that the convergence properties of the perturbation series are clearly 
poorer for K = 3 than ioT K = A suggests that there should be better solutions for the 
former case. If we try 

f,{x) = ^ c„2„+ia;2-+i, J = 0, 1, . . . (22) 

m=0 

then the values of log \Ej/Eq\ are smaller than those obtained earlier (compare K = 3{b) 
with K = 3 in Figure Q). The coefficient with the smallest absolute value is -E15 = 
-0.2344066313 x 10"^ and our best estimate resuhs to be El^^l = 1.1448015. 

It is well known that Fade approximants give considerably better results than power 
series J7| . We have tried diagonal Fade approximants [A^, N] on the perturbation series for 
the cases K = 2, K = 3 (b), and K = 4 and show results in Table HI Notice that the 
Fade approximants sum the K = 2 series to a great accuracy but they are less efficient 
for K = 3{b) and K = 4. This is exactly what is known to happen with the standard 
perturbation series for anharmonic oscillators 8]. However, Fade approximants appear to 
improve the accuracy of present perturbation results in all the cases discussed above. 

IV. CONCLUSIONS 

Fresent numerical investigation on the perturbation method proposed by Bessis and 
Bessis jl| suggests that although the series may be divergent they are much more accurate 
than those derived from standard Rayleigh-Schrodinger perturbation theory. One obtains 
reasonable eigenvalues for difficult anharmonic problems of the form V{x) = x"^^ , and re- 
sults deteriorate much less dramatically than those from the standard Rayleigh-Schrodinger 
perturbation series as the anharmonicity exponent K increases. In order to facilitate the 
calculation of perturbation corrections of sufficiently large order we restricted our analysis to 
polynomial solutions that are suitable for the ground state. The treatment of rational solu- 
tions for excited states (like those considered by BB [l|) is straightfordward but increasingly 
more demanding. 

Following BB [jj we have implemented perturbation theory by transformation of the 
linear Schrodinger equation into a nonlinear Riccati one. In this way the appropriate form 
of each potential coefficient Wj reveals itself more clearly as shown in Section IIIII for the 



quartic model. However, in principle one can resort to any convenient algorithm because 
the perturbation method is sufficiently general as shown in Section 2. 

A remarkable advantage of the method of BB [l|, which may not be so clear in their 
paper, is its extraordinary flexibility as shown by the two solutions obtained above for the 
case K = 3. Moreover, the method of BB, unlike two other renormalization approaches 
derived above as particular cases |3|,|j|, does not require and adjustable parameter to give 
acceptable results. 

We believe that further investigation on the perturbation method of BB [1] will produce 
more unexpected surprises. 
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FIG. 1: log \Ej/Eq\ vs. j for the ground states of the quartic {K = 2), sextic {K = 3 and K = "i 
(b)) and octic [K = 4) oscillators 
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TABLE I: Pade approximants [N, N] for the ground states of the quartic {K = 2) sextic {K = 3 (6)) 
and octic (K = 4) oscillators 



N 


K = 2 


K = 3{b) 


K = 4 


1 


1.06099633526211 


1.13079274107556 


1.210520271 


2 


1.06035870794451 


1.14501385704172 


1.203379654 


3 


1.06036222079274 


1.14479406990471 


1.225903753 


4 


1.06036204322721 


1.14480776322951 


1.225811686 


5 


1.06036210029130 


1.14480229840545 


1.225826780 


6 


1.06036209058862 


1.14480243489611 


1.225839331 


7 


1.06036209060731 


1.14480245393688 




8 


1.06036209048882 


1.14480245334362 




9 


1.06036209048295 






10 


1.06036209048246 






11 


1.06036209048301 






12 


1.06036209048423 






13 


1.06036209048420 






14 


1.06036209048417 






15 


1.06036209048417 






16 


1.06036209048408 






17 


1.06036209048418 






18 


1.06036209048418 






19 


1.06036209048418 






20 


1.06036209048417 






21 


1.06036209048418 






22 


1.06036209048418 






23 


1.06036209048418 






Exact 


1.060362090484183 


1.14480245380 


1.22582011 
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